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Asymptotic expansion for a solution of an elliptic boundary- value problem in a 

thin cascade domain U 



A. V. Klevtsovskiy, T. A. Mel'ny 



Abstract 

Asymptotic expansion is constructed and justified for the solution to a nonuniform 
Neumann boundary- value problem for the Poisson equation with the right-hand side 
that depends both on longitudinal and transversal variables in a thin cascade domain. 
Asymptotic energetic and uniform pointwise estimates for the difference between the 
solution of the initial problem and the solution of the corresponding limiting problem 
qq . are proved. 

<N . 
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S ■ 1 Introduction 



There are many articles and books (see, e.g., [T]-[T2]) devoted to boundary- value problems in 
thin domains (one of linear sizes of such a domain is substantially smaller than the others). 



\Q ■ The reason for such popularity of these problems is the wide possibilities of application of 

results in applied problems. In spite of a huge progress of computational tools it is impossible 
^ | to find acceptable numerical solutions of boundary- value problems in such areas because very 

■^ | small thickness of a thin domain naturally leads to a lengthening of computation time and 

significantly complicates the maintenance of an acceptable level of accuracy. Thus the main 
method of such researches is asymptotic analysis. The aim of this analysis is to develop 
rigorous asymptotic methods for boundary-value problems in thin domains. 

In recent years, the development of modern technologies of production of porous, com- 
posite, and other microinhomogeneous materials and biological structures has stimulated 
a significant interest in the investigation of boundary-value problems in thin domains of 
more complex structures: thin perforated domains with rapidly varying thickness and dif- 
ferent limit dimensions [131 E3], thin perforated domains with rapidly varying thickness 
[HI 031 H51 EE1 OH HE|, junctions of thin domains P3H [201 I2D 1221 123| , thick junctions of thin 
domains [2H 1221 1251 121], periodic grids and frames [2E1 I2§] - 

Research of various physical and biological processes in channels are topical for many 
branches of science (see, e.g., [301 [31] an d references indicated there). Of great interest to 
researchers is the appearance of different effects (such as sticking to walls, stenosis) in close 
vicinity to local irregularities (widening or narrowing) of geometry of channels. In [30l El] 
the author has summarized results of recent theoretical, experimental and numerical studies 
of flows and wall pressure fluctuations in channels with different types of narrowing. 
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Figure 1: Thin cascade domain 



In [23] by methods of formal asymptotic analysis, the limiting problem was obtained for 
a homogeneous Neumann problem for the Poisson equation with the right-hand side, which 
depends only of one longitudinal variable, in junctions of thin domains. The authors showed 
that the local geometric irregularity in the joint zone does not affect the view of the limiting 
problem. But, the convergence theorem and asymptotic estimates have not been proven. 

It should be stressed that the error estimate and convergence rate are very important 
both in justifying the adequacy of one-dimensional (two-dimensional) models to real thin 
three-dimensional models, and in the study of boundary effects and effects of local (internal) 
inhomogeneities in applied mechanics. Such estimates can be proved by developing of new 
asymptotic methods. 

In this paper we begin to develop asymptotic methods for boundary-value problems in 
thin cascade domains, which are cascade connections of thin domains with different thickness 
(see Fig. 1). To construct the formal asymptotic expansion we have generalized asymptotic 
method for thin domains with constant thickness proposed in monograph [12] . In particular, 
we introduced additional inner boundary-layer asymptotic expansion in a neighborhood of 
the joint of thin domains and studied its properties. Thus, the asymptotics for the solution 
consists of three parts: regular part, boundary part near the extreme vertical sides and inner 
part in in a neighborhood of the joint zone. 

It is understood that there is no principal difference between the construction of the 
asymptotics for the solution to a boundary- value problem in a thin cascade domain consisting 
of two thin domains of varying thickness, and in a thin cascade domain consisting of n thin 
domains of different thickness. Therefore, we consider a model two-stage thin cascade in this 
paper. Also, for simplicity, we consider two-dimensional case. 

The aim of this paper is to construct and justify the asymptotic expansion for the solution 
to a nonuniform Neumann boundary-value problem for the Poisson equation with the right- 
hand side that depends both on longitudinal and transversal variables in a thin cascade 
domain, which consist of two thin rectangles of different thicknesses e h± and e h,2 ■ 

The paper is organized as follows. In Section 2 we construct formal asymptotic expansion 
for the solution to the problem (TjQ). Section 3 is devoted to the justification of the asymptotics 
(theorem 13.11) and to the proving of the asymptotic estimates for the main terms of the 
asymptotics (corollary 13. ip . In the fourth section we analyze the results obtained and show 
possible generalizations. 



1.1 Statement of the problem 

A model thin cascade domain Q £ consists of two thin rectangles 

Of'^f-LOlxTf) and nJO = ((0, 1) X Tf ) , 



where Ye = (— £j,£j) , z = 1, 2; e is a small parameter; h± and /i 2 are fixed positive 
constants, h 2 < hi (see Fig. 2). 
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Figure 2: Model thin cascade domain fi £ 
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In £l £ — ((— 1, 0) xT £ ) (J ([0, 1) xT e 1 we consider the following mixed boundary- value 



problem 



-Au £ {x,y) = f(x, |), (x,2/)6 E) 

-9 v u e (x,j/)| ,,,3 = e<P±(x), x e Ii, i = 1,2, 

u e (-l,y) = 0, j/GTF, 

u e (l,3/) = 0, 2/GT £ 2) , 



[«e]U=0 
[<9x«e]U=0 



0. 



0. 



(1 



ye? 

yer 



yeri 2 , 



(2 



\T 



(2) 

£ 5 



where ii = (-1,0), i 2 = (0, 1) , [w]| x=0 = M(x,y)| x=0+ - M(x,y)| x=0 _ is the jump of func- 
tion u, d x = d/dx, d xx = d 2 /dx 2 . We assume that functions / and {(p± } are smooth 
functions in their domains of definition. 

From the theory of linear boundary- value problems it follows that for every fixed e there 
exists a unique generalized solution u e to problem ([!]), i.e., u e G H 1 ^,,), its traces on the 
vertical sides of fi £ are equal to {u £ \ x= ±i = 0) and the following integral identity 



/ Vu £ ■ \7ip dxdy = / / ' ip dxdy ^f e 2_\ / f± 4> dx 

holds for any function ip G H 1 ^^ such that ip\ x =±i = 0. 

Remark 1.1. In right-hand side of identity |lj) the following brief record 



(2) 



T^/2 / ^± ^^ x := 



2 f i 2 f 

y / <^+ ip dx + e 2_] / V 9 - ^ ^ 



zs introduced that will be used further. 

Our aim is to construct and justify asymptotic expansion of the solution u £ as e — > . 



2 Formal construction of asymptotic series 

2.1 Regular part of asymptotics 

Regular part of the asymptotics is sought in the form 

+00 

^^=Y. £k (^ )(x ^ + £ ~ 2 ^ )( ^)> (^y)enP, 2 = 1,2. 

fc=2 



(3) 



Formally substituting the series ((3]) into the differential equation and into the first boundary 
condition of problem (TJQ), we obtain: 



J2e k diut\ x , V ) - xy-'aj^M) - E 



d 2 u {i) v 

^ *■(*)« /(*,»*), V = y -, 



k=2 



k=2 



k=2 



dx 2 



+00 



E^ 






h; . 



2 (i) f \ 
£ tp y ±>(x) 



k=2 



Equating coefficients at the same degrees of e , we derive recurrence relations to determine 



,(<) 



coefficients {u k }. Let us consider the problem for u 



(0 . 



^2 „,(«) 



2, ,w 



tt^:')) = f( x ,V) + 



d 2 co 



.(0 



(*), 



<ix 2 



x 



r^GT, 



^< ; (x,r^)| .^ = </?±(x), xgZ 



(4) 



.(0 



«2 1 (^ •))' 



0, x G /,. 



Here Tj = ( — 2 ' 2 ) ' ( M ( x > '))xi := Jr u ( x ^ r l)df]^ i = 1,2. For each number i this is a 
Neumann problem for the ordinary differential equation with respect to the variable 77 G Tj ; 
variable x is considered here as a parameter. The last relation is added for the uniqueness of 
the solution. The solvability condition for the problem (j3j) gives us the differential equation 



for function uj. 



(0 . 



d 2 u^ 
dx z 



,(*)/ 



,W 



f(x,T])dr]-ifY(x) + (p_(x), x e Ii (i = l,2). 



(5) 



T, 



Let U2 be a solution of the differential equation (jSJ) (boundary conditions for this differential 
equation will be determined later). Then the solution of problem (jlj) is uniquely defined. 
To determine the coefficients u% , i = 1,2, we obtain the following problems 



-dLuffai) 



2, ,w 



<i 2 u; 



x 



.(*)/ 



-ax ( x ,v)\ v=± ^ 



dx 2 
0, x G Jj 

0, x G /,-. 



r?GT 4 



(6) 



Repeating previous arguments, we find that u^ = and 



2, ,0 



era;. 



e?x 2 



x G /,-, i = 1,2. 



(i) 

Let us consider boundary- value problems for the functions u k , k > 4, i — 1, 2 : 



~3j r, u k\XiV) 



.w 






0, a; G Ii 

0, 16 /,. 



(7) 



Assume that we determined all coefficients u 2 , . . . , 4-i> w 2 > • • • ■> u k-i °f the expansion (J3]). 
Then from the solvability condition for problem ((7|) it follows that 



dx 2 



;r 



d 2 x u k _ 2 (x,r])dr] = -d 2 x 



T, 



U 



w 

fc-2 



(x,n)dn = 0, 



i.e., Co^ is a linear function and 



W 2 / ,W 

— -§-(s) = 0, x G J*. 
OX 



(8) 



Remark 2.1. Ti/Z now boundary conditions for differential equations (TJP and |2P are wi- 
known. They will be found in the following stages of construction of the asymptotics. 

Thus the solution to problem (J7|) is uniquely determined, and hence the recurrent proce- 
dure for determining the coefficients of series fl3]) is uniquely solved. 

Remark 2.2. From the recurrent procedure of problems ([?]) it is easy to check that n^p+i 
are identically equal to for odd indexes k = 2p + 1, p £ N, i = 1,2. 



2.2 Boundary asymptotics near vertical sides of domain Q £ 

In the previous section the regular asymptotic expansion, which takes into account the right- 
hand side of the differential equation of problem ([1]) and the boundary conditions on the 
horizontal sides of the thin cascade domain Q e , was constructed. In this section we will 
build the boundary parts of the asymptotics, which neutralize the residuals from the regular 
parts of the asymptotics both on the left side of Vt e and the right one of fle . 

In a neighborhood of the left vertical side of fie we seek the asymptotic expansion for 
the solution in the form 

r (i) fl + x y* 



n« : =I>'n«(I±£f) 

fc=0 V ' 



(9) 



Substituting (j5j) into (jTJ) and collecting coefficients with equal degrees of e , we obtain the 
following mixed boundary-value problems 



rW 



A«x j (e,")= o, 



r(i) 



W ; (£^)l 



i ^1 



(e,7 ? )6(0,+c»)xT 1 , 

ee(0,+oo), 



rW, 



->«/ 



n^(o,n)= $^(n), neTu 
ni 1} (e,^)^ o, e^+oo, ne'T 



(10) 



1- 



,(!) 



,« 



-,(!) 



W/ 



,(!) 



where $^ = -«£ 2 (-l), fc = 0,1, Sj^fo) = -«£'(-!, 17) - < 2 (-l), k > 2. Here 



* = *?, u = *. 



Using the method of separation of variables, we find the solution of problem (ITU1) at a 
fixed index & : 



ni 1) K,,) = E 

p=0 



<)e"^ 



cos 



2 P^\ , di) --&£?*■* 



h 



, . , ...: -.■' ' ■ ; , / (2p+ 1)7T 

77 ] + 6£, >e h i ? sm ( ^ J — — —77 



:n: 



where 



,(D 



fci 



2p7T 



ill 
2 



*? W cos ( ^77 J d V , 6« = - y $« (77) sin , ] 



(2 P +1)7T 



r? d*7, 



,W 



2 2 



From the fourth condition in (TlOl) it follows that coefficient a must be equal to 0. As 
a result, we arrive at boundary conditions for the functions {u) k _l 2 } ■ 



c4+ 2 (-i) = 0, fee N . 



)(2) 



(12) 



In a neighborhood of the right vertical side of £l e we seek the asymptotic expansion for 
the solution in the form 



fe=0 \ 6 e / 



e e 

fe=U 

To determine the coefficients {11^ }fceN we get the following boundary- value problems: 
' -Ar v I$\?,r,)= 0, (r^)e(0,+oo)xT 2 , 

ni 1) (0,r / )= $f(r7), i]6T 2) 



(13) 



(14) 



where $i 2) = -wg 3 (l), fc = 0, 1; S^fo) = -uj, 2) (M) -w£> 2 (l), fc > 2; £* = ^ 
Similarly we find the solution of the problem ( 1T4"|) at a fixed index fc : 



+00 



nf(r,^) = E 

p=0 



C21 -^E2Lf* 

ai >e h 2 cos . , 



2/977 



(2p+l)n 

'e h 2 






(2p+l)7T 

sm I 77 



(15) 



where 



^2 h 2 



''2 
2 



,( 2 ) 



1 1 ^f^)^ = i y 4 2) (i^)^ - 4? 2 (i) = -<&(!)■ 



^2 ^2 

2 2 



From the fourth condition in (TH1) it follows that coefficient a must be equal to 0. It is 
possible when 

41(1) = 0, keN . (16) 



W = n f„„ U — O^ _l_ 1 ^ ^ T*J tn.^^^c *(*) 



a 



Remark 2.3. 5mce w^ = for k = 2p + 1, p G N, functions $2»+i = 0, p G N. ^4s 

nf = o, n^so, p gn, z = i,2. 

Moreover, from representation < f77]j and < f73|) ii follow the following asymptotic relations 

nj^&Tj) = 0(exp(-£O) *» e -> +oo, 4 2) (r,r/) = 0(exp(-£O) « s f "► +°°- ( 17 ) 

Equalities (TT2~j) and (TTBT) set boundary conditions at points —1 and 1 for all functions 
{u k } and {u;jj. } , respectively. We find out conditions for these functions at point in 
the next section. 

2.3 Inner boundary part of the asymptotics 

Let us consider what happens with the regular parts of the asymptotics in the join zone of 
two thin domains £l £ and £l e . Formally substituting the regular parts UoJ and UoJ into 
the transmission conditions of problem (fl|) , we obtain the following relations: 

+oo +oo 

£«* («?>(0, f ) +£ - 2 4"(0)) « £e l («i 2) (0, f ) + e -\»«(0) 

k=2 ' ' k=2 



JV ( 9^>(0, V - ) + ^(0) U JV 4«f(0, f ) + *-*£-«>) ) . (19) 



fc=2 \ / fc=2 

Equating, for example, the corresponding coefficients with the equal degrees of e in f|T8|) . 
we get 

u;f } (0) = u { k \o), when fc = 2 and fc = 2p+ 1; 

4 2) (0, f ) - 4 1} (0, |) = wg> 2 (0) - J k %(0), when fc = 2p, p G N. (20) 

The left-hand side of (12"U|) is a known quantity that depends on the "rapid" variable | and 

not necessarily equal to 0. Thus i1 
such that equality (12"0~|) is satisfied. 



not necessarily equal to 0. Thus it is impossible to choose the constant uJ k l 2 (0) — ^1+2^) 



Therefore, we should introduce an additional inner asymptotic expansion in a neighbor- 
hood of the joint zone to remove residuals that depend on the "rapid" variable - in the 

transmission conditions of problem ([I]) at the join zone of two thin domains Vt £ and Vt £ . 
Inner expansion is sought in the form: 



+00 

fc=i 



(21) 



Passing to the coordinates £ 



V 



in a neighborhood of the joint zone and then 



forwarding the parameter e to 0, we obtain the following unbounded domain 

S= ((-oo,0)xT 1 )u([0,+oo)xT 2 ), 

which is a union of semi strips E^ = (—00, 0) x Ti and S^ = (0, +00) x Y 2 . 

Let us introduce the following notation for parts of the boundary of the domain H : 

• <9S|| = {0} x (Yi \ T2) is the vertical parts of the boundary <9S , 

• dE& is the horizontal parts of the boundary 9SW, i — 1, 2, 

• dZ = = &E& U <92L 2) . 

Substituting (T2~Tj) into ([I]), taking into account residues that leave the regular parts of 
the asymptotics on the vertical sides and on the joint zone, and equating corresponding 
coefficients by equal degrees of e , we get the following relations for the coefficients {N k }. 

For even numbers k = 2p, peN: 



-AiV 2p 
d v N 2p 
d^N 2p 

W2 P ]k=0 

[d^N 2p }k =0 






e 2p 



2p 



■2p 



m a, 
on <9S = 
on 9S||, 
on T 2 , 
on T 2 , 



(22) 



where 



e^fo) 



duj 



(i) 

2p+l 



dx 



(0), 77 e as,. 



*2 P (v)=u$!(0,r))-u$(0,r)), V^T 2 , 



®2 P (V) 



d, , (1) d, , {2) 



dx 



For odd numbers /c = 2p + 1, p e N : 

-AiV 2p+ i 
d v N 2p +i 

[N 2p+ i]k=o 
[ [^iV 2p+1 ]| ?=0 








(0), r/GT 2 



m a, 
on <9S = 



2p+ i on <9H[| 
on T 2 , 



(23) 



$ 



2p+l 



on T 2 



where 



M)i 



dbj. 



(1) 

2p+2, 



e 2p+1 (7?) = -d x u^(0,rj) - I ^(0), V G AS,,, 

^ , {1) /// , (2) 

.W/'n ^ _ a „.( 2 Vn ^ i 2 p +2 (0) - 2p+2 (Q), r/ G T 



$ 2p+ i(??) = d x u%(0, ri) - d x u%(0,ri) + 



2; 



dx dx 

It should be noted here that u = 0, iti = 0. 

In order to find out whether exist functions that satisfy the relations of problems (1221) 
and (1231) . at first we study the solvability the following boundary- value problem: 



-AN(£,r)) 


= Ffotj), 


(e^)eS; 




d v N(^ V )\ v=± ^ 


= ±si°(0, 


(-l)<£e(0,+oo), 


z = 1,2 


^jv(e,»?)ic=o 


= G( V ), 


ieTi\T 2 ; 




Mle=o 


= *(v), 


r/GT 2 ; 




Mk=o 


= <%), 


r?GT 2 . 





(24) 



Let Co^(5) be a space of infinitely differentiable functions in H that are finite with 
respect to £ , i.e., 

\/veC^(E) 3R>0 V((,i,)GH |£|>fl: u^t?) = 0. 



Define the following space H := (C^(S), || • ||%) , where 



MN 



Vv(C,v)\ 2 d^dr ] + / K£,r/)| 2 |p(£)| 2 c^, 



and the weight function p(£) = (1 + |^|) _1 , (el, 

The case \l/ = . Let us take an arbitrary function v G C^ (5) , multiply the differential 
equation of problem ( l24"j) by t>, and integrate this equality over domain S . Using the Green- 
Ostrogradsky formula, we deduce the following integral identity: 

//■ r0 r+oo 

VN-Vvdidr,= / FvdZdr l + / fl^) v(£, ±^) d£ + / ^(O v(£, ±f ) d£+ 
J J-oo J0 

+ / G(r])v(0,r])d7]- ^(7])v(0,r])dr]. (25) 

Ti\T 2 T 2 

Definition 2.1. A function N from the space % is called a weak solution of problem (24\ ) 
if the identity (25\) holds for all v G % . 



From lemma 4.1, remark 4.1 and 4.2 (see [25]) it follows the following proposition. 

Proposition 2.1. Let p~ x F G L 2 (S), p^B^ G L 2 (0,+oo), p^B^ G L 2 (-oo,0), G G 
L 2 (Ti\T 2 ) and $ G L 2 (T 2 ). 



Then there exist a weak solution of problem [24\ ) if and only if 

F d£d V + / B^ (0 d£ + / 4 2) (0 d£ + / G( V )d V . (26) 

J-oo JO J 

T 2 3 Ti\T a 

In addition, this solution is defined up to an additive constant and we can choose this 
constant such that there will exist a unique solution of the problem [£$ with the following 
differentiable asymptotics : 

{C(exp(^-0) when f -»■ — oo, 
1 VK ^ ,J (27) 

d + 0(exp(-££)) when £ -> +oo. 

Moreover, if the functions F, G, $ are even with respect of rj (F, G, $ are odd with 
respect of rj) and B = B_j! , i — 1, 2 (B_ = —B+, i = 1,2), i/ien solution Ao zs even 
(ogMJ function with respect of r\. If Ao «s odd function, then in (2l\ ) the constant do is equal 
to zero. 

From corollary 4.1 ( [25J ) it follows the second proposition. 

Proposition 2.2. There exists a nontrivial solution Z of the corresponding homogeneous 
problem (jE$\ ), which does not belong to the space "H, and this solution has the following 
differentiable asymptotics: 

( |- + C hl + 0(exp(£O) when Z^-oo, 
Zo(Z,v)={ J 11 (28) 

j- + C h2 + 0(exp(-£O) w/ien ^ -»> +oo, 

w/iere C^ = /^(^(O,^))^, z = 1,2. 

Moreover, the function Zq is even with respect of variable r\ and any other solution of 
the homogeneous problem ( fl%y , which has polynomial growth when £ — > ±oo, zs a linear 
combination a\ + a^^o, where a.\ and a.2 are some constants. 

Remark 2.4. Using the second Green- Ostrogradsky formula, similarly as was done in re- 
mark 4-3 (JMEj), constant d from pT7| ) can be found as follows 

/pO />+oo 

F{z, v )Zo(z,v)dtdTi+ / 4 1) (0^o(^,±^)^+ / 4 2) (0^o(e,±f)^+ 
</-oo JO 

+ y G( V )Z (0, V )d V - J$( V )Z (0, V )d V . (29) 

Ti\T 2 T 2 

In the general case when \l/ 7^ 0, the following substitution 

W = N - X 5^ 
must be done in the problem (124")) . where xs £ C°°(K1+)) < ^ < 1 and 

{0, when £ < 0, 
1, when <£<<$, (<J>0) 
0, when £ > 25. 

Then [W]|f=o = and we arrive to the previous case. 

10 



Definition 2.2. Let \1/ e H 1 (T2). A function N is called a weak solution of problem [24\ ). 
if there exists a function W from 7-L such that the following integral identity 



jvW-Vvdidr ] = j Fvd^- j V{xs(^(r l ))-VvdZdr l + j^ b£ 1} (0 v(£,±^) <£+ 

B£ 2) (0«(£,±^)#+ J G(rj)v(0,ri)dTi- J $( V )v(O lV )dr) (30) 

Ti\T a T 2 

holds for all v £% . 

Remark 2.5. It is easy to verify that in the general case, the proposition \2. 1\ holds with the 
same solvability condition /[2b]) . Moreover, the function ^ must be even (odd) in the last 
item of Proposition \2.1\ 

Now we back to problems ( )22l) and ( 1231) . From (126]) it follows that the following equalities 

j e k {r])dr] = [ Mv)dv, keN, (31) 

03|| T 2 

are the corresponding solvability conditions for those problems. 

Taking into account that (u k (x, •))y i = 0, x G /,, i = 1,2, from ( 15T]) we derive the 
following relations for functions {ojj. } : 

A, .(1) J. .(2) 

/H^-(O) = /i2^-(0), keN, k>2. (32) 

Hence, if {u k %> } satisfy ( 1521 . then there exist solutions of problems ( 122]) and (125]) . According 
to Proposition 12.11 we can uniquely choose those solutions such that they have the following 
asymptotics 

\ d+ + e>(exp(-§0) when 5 -> +00. 
In what follows, in asymptotic expansion (J21I) we will use the functions 

N k (£ ,ri) = < keN. 

Then taking into account (155j) . functions {iV fc } are exponentially decrease as £ — > ±00. 

Formally substituting the series ([5]) and (12"T|) in the first transmission condition of problem 
dTJ, we get 

(«g + iVoo) | x=0 _ = ( w g + iVoo) U =0+ , 

from where, by virtue of the corresponding equalities in problems ( )22l and ( |23l . we deduce 
the following relations for the functions {u% } : 

w? ) (0)=w? ) (0), wi 2) (0)-wi 1) (0)=4_ 2 , fc6N, fc>3. 

11 



Thus, we have obtained a sequence of boundary-value problems to determine functions 
{uo k }. For functions u 2 and co 2 that form the main term of the regular asymptotic 
expansion fl3]), the problem looks as follows 



, d 2 uf 
-hi 



x) 



hi 



dx 2 

^(O) 

^ 1} (o) 



F«(x). 

4 2) (o), 



X G ij, 2 



1,2, 



dx 



ft. 



c/w. 



(2) 



d;r 



-(0), 



0, u4 2) (l) 



where 



p(0, 



•n 



,W 



(0/ 



f(x,7])dr]- (pY(x) + (f_(x), x e Ii, i = 1, 2 



The problem (1341) will be called homogenized problem for problem ([T]). 
For next functions {u k , co k . k > 3} we get the following problems: 



x 



h f«f 
dx 2 

«f(0) 



w 



c/.r 



0. 

4 2) (o)-4- 

o, 4 2> (i) 



x G Ij, i 



1,2, 



2' 



0. 



It is easy to find the unique solution to problem ( 136]) at a fixed index k : 



4 1! w 



fr-2 d k _ 2 
hi + /i 2 



'x + 1), x G Ii; w 



(2), 






1-x), x e l 2 - 



(34) 



(35) 



(36) 



(37) 



3 Scheme of construction of the complete asymptotics 
and its justification 



Let us introduce the following notation 



u k {x,-) 



u 



(i) 



u 



(2) 



(x, §), fay) e SI 

(x, f), (i,i/)6fl 



(i) 



X 1 ) 



(2) 



w fc (x) 



u k ; (x), x G A 



,(2) 



fceN, k > 2. 



tof. '(x), x G I 2 



From homogenized problem (1341) we uniquely determine the main term of the asymptotics 
0J2 of the series ([3]). Then from problems (j4]) that can now be rewritten as 



-dl v uf{x,r,) = f(x,rj)-h i 1 F^(x), r] G T u 



.«, 



-d v uy(x,r})\ v=± h± 



,W 



V9±(x), xG/j 

0, X G /,;, 



z=l,2, 



(3? 
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we uniquely determine 

v $ ) (x,rj) = - f (r 1 -t)(f{x,t)-hT 1 F {i) {x)\dt- W {i) {x)+af{x) 

where function a 2 are uniquely determined from third condition in ([38]) . i.e. 



(39) 



OCo (x) 



T, 



(r ] -t)f{x,t)dtdri-6- 1 h 2 i F ii \x), i = 1,2; 



functions F^ and F^ are given by formulas (135]) . 

Now with the help of formulas (TTTj) and (TT5l) , we find the first terms IT;, and fl 2 of the 
boundary-asymptotic expansions © and (IT3"|) respectively; they are solutions of problems 
(110 1) and (1141) that can be rewritten as follows 



-A^ 1) (e,r / )= 0, (£,77)e(0,+oo)xTi, 

An? ) (e^)l fl _±ii= o, £g(o,+oo), 

^'(0,1])= -^(-l^), ^T, 

n?^)-* o, e^+oo, fieTi, 



(40) 



r(2) 



-A e ,nr j (r,r/)= o, 



r(2) 



-d v U^(e, V )\ 



v =±- 



0. 



(e,^)G(0,+oo)xT 2 , 

ee(o,+oo), 



n( 1) (0,r / )= -t4 2) (1,77), r/GT 2 , 
n 2 1} (r,r7)^ 0, f->+oo, ??eT 2 . 

Then we find the first term of the inner asymptotic expansion (J2"T|) 

M (£,??), £<0, 

Mi{£,v)-dti e>o, ' 

where M is the unique solution of the problem ( 12"3"]) that can now be rewritten as 

-AM = in 2, 

dffhfx = on 3EL=, 



(41) 



#i(£,*7) 



<%M 

[M]| ?=0 

[%A/i]k=o 



dug 

dx 



(0) 







on <9S||, 
on T 2 , 



(42) 



-(O)--tMO) onT 2 , 



dx 



c/.r 



with asymptotics (133]) . Recall that constant g?|" is also uniquely determined (see remark [2~4]) . 
Thus we have uniquely determined the first terms of the asymptotic expansions ([3]), ©, 
(TT3D and (J2II). 
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Assume that we have determined coefficients cu 2 , ■ ■ ■ , Co>2n-2, «2, «4, • • • , W2n-2 of the series 
([3]), coefficients 11^ . , II4 , . . . , U. 2 l_ 2 of the series (Q and ([TBI respectively and coefficients 
N 1 ,..., N 2n _ 3 of the series (I2TJ) . 

Then, using formulas (13"T|) . we write the solution c^n-i of problem f[3T)j) with the constant 
d 2n _ 3 in the first transmission condition. Further we find the coefficient 



N 2 n-2^,V) 



Af 2n -2(£,V)-dt n -2, £ > 0, 



of the inner asymptotic expansion ( 12T]) . where N2n-2 is the unique solution of the problem 
fl23l) that can now be rewritten as 



-AA/2„-2 
d v Jv2n-2 

9^2n-2 

[A/2n-2J|f=0 

[<%A/2n-2]|£=0 






^2 ^2n-3 



m c, 
on r3E = , 

on9S||, 
Ail + ft 2 

4n- 2 (0, 77) - ^£-2(0, r?) on T 2 , 
dtn-aihi - h 2 ) 



(43) 



/ii + h 2 



on T 5 



and J\f 2n - 2 has the asymptotics ( 1331 . 

Knowing G?2 n _ 2 and using fl37|) . we get the solution u 2n of problem f )36l) . Next coefficient 



N^iforj) 






of the inner asymptotic expansion ( |2J1 is defined with the help of solution ftf 2n -i to problem 
(1231) that can now be rewritten as follows 



<%A/2n-l 

[AT 2 n-l]k=0 

[<9 5 A/2n-l]k=0 








-fli B «S- a (o^) + / ' 24 "- 2 







,« 



/ll + /l2 



(2) 



<Wn- 2 (0, 77) - <9x«2n_ 2 (0, //) + 



d tn-2( h l ~h 2 ) 



in a, 

on <9S = 

on 9H|| 
on T 2 , 

on To. 



(44) 



Coefficients m 



2n; 



/ii + h 2 

= 1,2, are determined as solutions of following problems 
-9^ v u 2n( x ,v) = d% x u { 2 i_ 2 {x,ri), rjeTi, 

-dr,U ( 2 i(x,T))\ v=± h± 
(«£(*, -)>T, 



0. 
0. 



x e Ii 
x e h. 



(45) 
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And finally, we find the coefficients U. 2 n and EE^ of the boundary asymptotic expansions 
(|HJ) and ( Tl3l) respectively as the solutions of problems (TTU1) and (tHJ) that can be rewritten 
now as follows 

-A^n2(e,r/)= 0, (^»|)€(0,+oo)xTi, 

-^n£(e^)L ±il = o, eG ( ,+oo), 

£ -> +oo, »7 e Ti, 



■«o 



n£ ) (o,r ? ) = 

' -A e?? ng(e,r/)= 0, (e,^)G(0,+oo)xT 2 , 

-a t| n£ ) (r,^)l t|=± ^= o, fe(o,+4 

nS(0,t/)= -tig(M), r/eT 2 , 
nS(e*,^)^ 0, T^+oo, r/GT 2 . 

Thus, we can successively determine all coefficients of series (J3J), ©, (|13l) and (T2TT) . 



(46) 



(47) 



With the help of the series 



(fT5|) and (JUJ) we construct the following series 



w 2 (x 



+oo / \ +oo / 

+ J2 ^ " 2k+l {x) + X °(x)N 2k ^ (-, -) + X) ^ M2fc ( x ' ") + ^ +2 (x) + 

*;=i V / fe=i V 



(i)/l + x y\ +( {2 )/l-x y 



+ A*)N 2t (f,f) + x»n« (i±£, f ) + x+( ,)nS> ( 



(i,l)Ell„ (48) 



where x , X are smooth cut-off functions defined by formulas 



X ± {x) 



1, if |1 + x| < 8, 
0, if |l + x| > 25, 



X°(x) 



1, if |x| < 8, 
0, if Ixl > 25. 



Here 5 is an arbitrary sufficiently small fixed positive number. 



Theorem 3.1. Series ( j 48 1) zs i/ie asymptotic expansion for the solution of the boundary- 
value problem ( IT]) m £/ie Sobolev space H l (Q £ ). Moreover, the following asymptotic esti- 
mates 



VmeN 3C m >0 3e >0 \/ee{0,e ): 
hold, where 

m 

Ui m \x,y) = u 2 (x)+J2z 2k ~ l (w 2 fc+i(x) + X°(x)N, 



U f 



-U^\\ H , m <C m e 2m ^ (49) 



2fc-l 



fe=l 



e e 



in I 

+ ^2 £2k [ U 2k [X,-J+U 2k+2 (x) 



+ ^(*)^(p|)+x>^ (s.iOen., (so) 



is the partial sum of ( 1 48 



15 



Remark 3.1. Hereinafter, all constants in inequalities are independent of the parameter e. 

Proof. Consider an arbitrary m£N. Substituting the partial sum Ug into equations and 
boundary conditions of the problem (TjQ) and considering relations (|34l)-fH71) that satisfied by 
the coefficients of the series (HHD, we find 



AU^ m \x, y) + f (x, |) = e 2m d xx u 2m (x, |) + 

171 J j2 

+ ^, 2fc - 1 (2,-^(x)^iV 2fc _ 1 (^r / ) + ^(a;)iV 2fc _ 1 (^r / ))| 5=f , r;=f + 



m A 

+ ^2e 2k {2e- 1 ^(x)d i N 2k ^ V ) 

k=l 



2,_0 



dx 2 



{x)N 2k (C,v) + 



Taking into account the exponential decreasing of functions {iV^, 11^. , 11^. } (see (|33|) 
171)). we conclude that 



(51) 



3 C m 3 e \/ee(0,e ): lim sup | R™ (x,y)\< C m e 2m (i = 1 , 2) . (52) 

It is easy to check that the partial sum leaves the following residuals 



r( m )/ 



,W/ 



-d v U} '(x,y)\ y=±£ ^ = eip ( ±>(x), 
Ut\±l,y) = 0, 

d x ut\x,y)\ x=0 = e 2 ™ (d xU U (0,*) - j^- d: 

in the boundary conditions (i — 1, 2) and the following ones 



2m 



y g if, 



=:*T(y), »6W 



£ 5 



[tf« 



Mil 

£ |x=0 



0. 



1/ e TQ ; , 



r( m )l 



h 2 - hi 



>M/ 



,(2) 



[^ni,= ^([^^u+^q^o^^r^y), i/gt 



in transmission conditions. Obviously that there exist positive constants C m and Eq such 
that 



\/ee(0,e ): lim sup i^ m) (y) 

^ eT (l)^ T (2) 



< C m e , lim sup 



yer 



(2) 



^ m) (y) 



S C m £ 



2m 



(53) 
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Thus, the difference W, 



W £ :=u e - m 


satisfies the following system: 


-aw £ 


_ fi( m ) 


in Q £ , 


d y W £ (x,±e^) 


= o, 


x e I t , i = 1,2, 


W e (±l,y) 


= o, 


yeT £ l \ z = l,2, 


-d x W e (0,y) 




y G T^Tf), 


[W e }\ x=0 


= o, 


r- -Y-( 2 ) 


[d x W £ }\ x=0 


p( m ) 
— —tie , 


ye TV. 



(54) 



This means that the series (148]) is a formal asymptotic solution of problem (flj) . 
From (15411 we derive the following integral relation: 



f \VW £ \ 2 dxdy = f R^ W £ dxdy - J R { ™ } (W £ \ x=0 ) dy + J R^ (W £ \ x=0 ) dy . 



<>- 



y( 1 )\ y^ 2 ) 



-(2) 



Now, using the Friedrichs inequality and estimates ( 152]) and ( 153]) . we deduce from the previous 
integral equality the following inequality 



/ \VW £ \ 2 dxdy < C m ^hi+h 2 e 2m+1 ? \\W E \\ L 2 m + C m \[h x e 2m+ ^ 



• 5 ||^ £ (0,-)IL 2(T ( 1 ) ) < 



<C m e 2m ^\\VW £ \\ LHn£) , 



which means the asymptotic estimates (|49|) are satisfied. Those asymptotic estimates justify 
the constructed asymptotics and imply that series ( 1 48 j) is the asymptotic expansion for the 
solution of problem (jTJ in iJ 1 (r2 £ ). □ 

Corollary 3.1. For difference between solution u £ of problem (QJ) and solution u 2 of the 
homogenized problem [34\ ) the following asymptotic estimates 



\u £ -u) 2 \\tf((i e ) <C ei, \\u £ -u 2 \\m(n E ) < C e, 
\\u £ - w 2 - e{u 3 + X°Ni)\\ Hl{nE) < Ci e'i, 



\\E^(u £ )-u 2 \\ LHu) <C 2 e, 



\Ei l \u £ )-u 2 \\ m(h) <C 2 eK 2 = 1,2, 



hold, where 



max 
xeii 



\Ef{u £ ){x)-uj 2 {x)\ <C 3 e*, 2 = 1,2, 



Ef\u £ ){x) 



e hi ly(i) 



u £ (x,y)dy, 2 = 1,2. 



(55) 

(56) 

(57) 
(58) 



i JT 
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Proof. Since functions iV 2 , n^ and IT 2 exponentially decrease at infinity, from estimate 
( )49l at m = 1 we deduce inequality ( )56l) : 

||m c - w 2 - e(w 3 + X°Ni) || H i (ne) < 

< Ik - 0. (1) llfli ( a.) + IN 1} - W2 - ^ + x M)L 1(Oe) < 



< Ci £2 + £ S 



M 2 + ^4 + x°iv 2 + x"n^ ) + x + n? 



r(2) 



H" 1 (n s 



< 



<£ z 



v4n 2 + u A + x °iv 2 + x-r4 x) + x + nf) 



+£' 



W2 + u; 4 + /iv 2 H- x -n( 1) H-x + nf 



L 2 (Q E ) 



+ 
£ 2 (fie) 

+ C 5 ! el < 



t=i 



< e 2 ^ ( / |^u 2 (x,?])| 2 (ix(ir7J 2 + 



/,xT, 



2 
+^ 5 5Z ( C 2 1 1^4 1 1 ^1(7,) + ||w 2 |Ua(j. xT .) + ||9 ;r M2||L2( /iXTi) J + 

+ e ( IIv^a^Hl 2 ^) + l|v^n 2 ||x 2 ((o,+cx3)xY 1 ) + l|v^*^n 2 ||l 2 ((o,+oo)xt 2 )J + 



+ £ 3 (||^V 2 || L 2 (H) + ||n^ 1) || L 2 ((0i+oo)xTl) + ||n^ 2) || L 2 ((0i+oo)xT2) j +de'i <Ci£ 2 . (59) 

Taking into account that ||x°-Wi||r, 8 (fl e ) < c 2 e and llx^ill^imW-, < c 3 5 i = 1,2, from 
(15E|) we get (155]) . 

Using the Cauchy-Buniakovsky inequality, from (155]) we derive inequalities fl57|) . Since 
the space i? 1 (/i) continuously embedded in C(Ii), from the second inequality in (j3"7j) it 
follows inequality (158]) . D 



- 3 

Remark 3.2. Summands of order e^ in the sixth line of (I 59 \ bring the main contribution 
to the determination of the constant C\ in inequality (I 56 j) . Taking the explicit form of the 
coefficients u 2 {see (I39p ) into account, we can specify the dependence of this constant on 
the right-hand sides of problem ( IT"}) : 



y2( / \d v u 2 {x,r])\ 2 dxdi]) 2 < 

2 



i=l 



-ill/Hi^/iXTi) 



+ 2^11^11^) H-V6||^ || La(/4) ). (60) 



Terms of order e 2 in the eighth line of (1 59 1) import next contribution to the constant 
C\. Let us estimate those terms. From the corresponding integral identity for the solution 
N 2 (see ( j30p ). we derive 

l|V^JV 2 || L2(a) < c(5)||* 2 || L2(T2) + \\%\\LHr 2 ) + 2di h2 ^f+ hiy/r2 < 

ill ~r tl2 
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< 2J2 ( ^(c(5)v / 2^ V / l7^+ VS) (||/||z»(7,xT,) + \\d x f\\ L 2 {IiXri )) + 



+ y/h i (c{5)h i ylz + 2h 2 i +y/£) (\\<p® || L2(/i) + \\d x <$ || L2(/i) ) ) + 

+ ^ fe 1 + fe 2 ' (61) 

where ^O?) = m 2 (0,n)— it 2 (0, 77). Similarly we can estimate values HV^ID, ||l 2 ((o,+oo)xTi) 
and ||V^II 2 ||i3((o,+oo)xT 2 )- 

Remark 3.3. // p± = and function f in the right-hand side of the differential equation 
of problem (jTJ depends only on the variable x, then all coefficients {n 2 fc}, {n 2fc } and 



r(2) 



{II 2fc j } are equal to 0. In this case the asymptotic series ( 1 48 j) /10s the following fi 



orm 



+00 

fc=l 



(x) + Y, £k [ "W*) + X°(x)N k (p I) , (x, j/) e fi e , (62) 



and asymptotic estimate (I 56 1) zs of order e 2 . Moreover, the value HV^A^H^^), which is 

{hr±hy\ 

hj+h-2 



now bounded by the value 2a'j~ ..^IV (see (1 6 1 j) ) . frrmg i/ie mam contribution to the 



constant C\ in f l 56 j) . 

4 Discussion of the results 

1. Estimate (I56"j) shows the structure of the corrector in the asymptotic approximation for 
the solution u £ of the problem (GQ). It (corrector) has the form 

eM*) + X°(aOtfi(f.f)). 

and its gradient is equal 

X°(x)V^iV 1 (^r / )| 5=f ,, =f +, (^(x) + (x°(aO)'tfi(f, f)) . 
Since function N\ exponentially decreases at infinity (see (j3"3"|) ). then 

e ll w 3 + X ^i|| H i (n , ) <Coe- 

Thus, for this structure of thin cascade junctions there are not any significant boundary 
effects in a neighborhood the join zone for the solution u £ of problem (Tj[|). This also indicates 
second estimate in ( 155]) and uniform pointwise estimate fl58|) . 

The results obtained give the right, in terms of practical application, to replace the 
complex boundary-value problem ([T]) with the corresponding simple boundary- value problem 
with sufficient accuracy that measured by the parameter e characterizing the thickness. 
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Moreover, in this paper we show how the constant C\ in the main asymptotic estimate 
( |56|) depends on the right-hand sides of problem ([1]) and on the geometrical parameters h\ , 
]\2 and df (see remark I3.2J) . Also it is possible to indicate how other constants in the 
asymptotic estimates from Corollary 13 . 1 1 depend on these values. This fact makes it possible 
to directly use the asymptotic estimates for the approximation of solutions of boundary- value 
problems in thin cascade domains instead of numerical calculations. 

2. The method proposed in this paper for the construction of asymptotic expansions can 
be applied without substantial changes to the asymptotic study of boundary value-problems 
in thin cascade domains with more complex structures, namely, either thin cascade domains 



£/?, 



Eh. 



V 



0(8) 



Figure 3: Thin cascade domain with local widening 



with local widening (narrowing) in a neighborhood of the join zone (see Fig. 3), or thin 




Figure 4: Thin cascade domain of graph type 



cascade domains of graph type (see Fig. 4), or thin cascade perforated domains with rapidly 



000000000000 



Figure 5: Thin cascade perforated domain with rapidly varying thickness 



varying thickness (Fig. 5). We need to add series with rapidly oscillating coefficients (see 
[T5| [T4] ) to the regular part of the asymptotic expansion for solutions of boundary- value 
problems in thin cascade perforated domains with rapidly varying thickness. 
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